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We present a measurement of time-dependent CP asymmetries and an updated determination of 
the CP-odd fraction in the decay B° — > D* + D*~ using a data sample of 88 x 10 6 BP> pairs collected 
by the BABAR detector at the PEP-II B Factory at SLAC. We determine the CP-odd fraction to be 
0.063 ± 0.055(stat) ± 0.009(syst). The time-dependent CP asymmetry parameters Im(A+) and |A+| 
are determined to be 0.05 ± 0.29(stat) ± O.lO(syst) and 0.75 ± 0.19(stat) ± 0.02(syst), respectively. 
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The Standard Model predicts these parameters to be — sin2/3 and 1, respectively, in the absence of 
penguin diagram contributions. 



PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 

The symmetry for combined charge conjugation (C) 
and parity (P) transformations is violated in B decays. 
Measurements of CP asymmetries by the BABAR Q and 
BELLE 0] collaborations established this effect and are 
compatible with the Standard Model expectation based 
on the current knowledge of the Cabibbo-Kobayashi- 
Maskawa Q quark-mixing matrix. As a result of the 
interference between direct B decay and decay after fla- 
vor change, a CP-violating asymmetry is expected in the 
time evolution of the decays B° — > D* + D*~ Q within 
the framework of the Standard Model || . This CP asym- 
metry is related to sin2/3 when corrections due to theo- 
retically uncertain penguin diagram contributions are ne- 
glected . Penguin-induced corrections are predicted 
to be small in models based on the factorization approxi- 
mation and heavy-quark symmetry; an effect of about 2% 
is predicted by Ref. Q- A comparison of measurements 
of sin2/? from b — > ccs modes such as B° -» J/^K% @ 
with that obtained in B° — ► D* + D*~ is an important 
test of these models and the Standard Model. 

The B Q — > D* + D*~ mode is a pseudoscalar decay to 
a vector-vector final state, with contributions from three 
partial waves with different CP parities: even for the S- 
and D-waves, odd for the P-wave. The CP-odd contri- 
bution is predicted to be about 5.5% in Ref. 0. We 
present an updated determination of the CP-odd 
fraction, R±, based on a one-dimensional time- integrated 
angular analysis. We also present a measurement of the 
time-dependent CP asymmetry, obtained from a com- 
bined analysis of the time dependence of flavor-tagged 
decays and the one-dimensional angular distribution of 
the decay products. The data used in this analysis were 
collected with the BABAR detector at the PEP-II stor- 
age ring. The BABAR detector is described in detail 
elsewhere [l^. The data sample corresponds to about 
88 x 10 6 e+e~ -> T(45) -> BB events. 

B° mesons are exclusively reconstructed by combin- 
ing two charged D* candidates reconstructed in the 
modes D*+ -> D°n+ and D*+ -> D+n°_ We in- 
clude the D* + D*~ combinations (D°ir + , D°ir~) and 
(D°7T+,p-7r ), but not (P+tt , B^tt ) due to the smaller 
branching fraction and larger backgrounds. Prior to 
forming a B°, the D* candidates are subjected to a mass- 
constrained fit and vertex fit that includes the position 
of the beam spot. 

The reconstructed D° and D + modes are D° — > 
K~ir + , K~ir + 7r°, K~ir + 7r + ir~ , K®tt + tt~, and D + — » 
AT~7r + 7r + , ATg7r + , K~ K + ir + . The reconstructed mass 
of the D° (D + ) candidates is required to be within 
20 MeV/c 2 of the nominal D° (£>+) mass H3, except 
for D° — » K~Tr + Tr Q , which has a looser requirement of 



35 MeV/c 2 . The D candidates are subjected to a mass- 
constrained fit prior to forming D* candidates. 

Charged kaon candidates are required to be incon- 
sistent with the pion hypothesis, as inferred from the 
Cherenkov angle measured by the Cherenkov detector 
and the specific ionization measured by the charged- 
particle tracking system. No particle identification re- 
quirements are made for the kaon from the decay D a — > 
K~ir + . The reconstructed mass of K® — > 7r + 7r~ candi- 
dates is required to be within 25 MeV/c 2 of the nominal 
K° mass. The angle between the flight direction and 
the momentum vector of the K® is required to be less 
than 200 mrad, and the transverse flight distance from 
the primary event vertex must be greater than 2 mm. 
A mass-constrained fit is applied to each K® candidate. 
Neutral pion candidates are formed from two photons de- 
tected in the electromagnetic calorimeter, each with en- 
ergy above 30 MeV; the mass of the pair must be within 
20 MeV/c 2 of the nominal ir° mass, and their summed en- 
ergy must be greater than 200 MeV. A mass-constrained 
fit is applied to these ir° candidates. The mass of the 
7T° from D* + — > P + 7T°, however, is required to be within 
35 MeV/c 2 of the nominal 7r° mass, and the momentum 
in the Y(AS) frame in the interval 70 < < 450 MeV/c, 
with no requirement on the photon energy sum. 

We construct a mass likelihood Cyiass that includes the 
mass and mass uncertainty of the D and D* candidates. 
The D mass resolution is modeled by a Gaussian whose 
variance is determined on a candidate-by-candidate ba- 
sis. The D* - D mass difference resolution is modeled 
by a double-Gaussian distribution whose parameters are 
determined from simulated events. The value of CMass is 
used to select B° candidates, with a different requirement 
used for each D decay mode combination. In an event 
where more than one B° candidate is reconstructed, the 
candidate with the largest Cmuss value is chosen. 

The primary variables used to distinguish signal from 
background are the energy-substituted mass, toes = 
V -^icam — Pb ' an d the difference of the B candidate 
energy from the beam energy, AB = Eq — BB eam , 
where all variables are evaluated in the T(45) center- 
of-mass frame. The B° candidates are required to have 
-39 < AB < 31 MeV and to es > 5.2 GeV/c 2 . 

To reject backgrounds from the e + e~ — > cc continuum 
process, events are required to have a ratio of second 
to zcroth Fox- Wolfram moments ^4| of less than 0.6. 
We also require that the cosine of the angle between the 
thrust axis of the reconstructed B and the thrust axis of 
the rest of the event be less than 0.9. 

After all selection criteria have been applied, a fit to 
the toes distribution using a Gaussian and an ARGUS 
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function [lfj for the signal and background, respectively, 
results in a signal yield of 156 ± 14(stat) events. In the 
region m ES > 5.27 GeV/c 2 , the signal purity is 73%. 

We perform a one-dimensional angular analysis to de- 
termine the fraction, R±, of the P-wave, CP-odd com- 
ponent of the B° — > D* + D*~ decay. In the transversity 
basis the following three angles are defined: the an- 
gle 0\ between the momentum of the slow pion from the 
D*~ in the D*~ rest frame and the direction of flight of 
the D*~ in the B rest frame; the polar angle #tr between 
the normal to the D*~ decay plane and the direction of 
flight of the slow pion from the D* + in the D* + rest 
frame; and the corresponding azimuthal angle </>t r . The 
time-dependent angular distribution of the decay prod- 
ucts is given in Ref. [lfil ]. 

The dependence of the detector efficiency on the de- 
cay angles can introduce a bias in the measured value of 
R±. Including the efficiency explicitly in the decay rate 
and then integrating over time and the angles 9\ and 4>tc 
results in the one-dimensional differential decay rate: 



1 dT 



T d cos 9 t 
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[(1 - R ± )sm 

- Ct 

— J o (cos0 tr ) 



2R i cos tr x 1 1 



I- a 

(cOS0 tr )] 



I\\ (cos tr ) 



(1) 



where R± 
M|)/(M 2 . 



M 2 /(M 2 



M \\ 



Ml) 



(M 2 



My ), and (M , M\\, M±) are the magnitudes 



of the amplitudes in the transversity basis. The three 
efficiency moments, Ik (k = 0, ||, _L), are defined as 



Ik (cos 9 tI ] 



efcos 6»i d^tr 9h (0i , 0tr ) e {9i , 9 tr , (f>tr ) , (2) 



where go = 4 cos 2 9\ cos 2 </>t r , g\\ = 2 sin B\ sin tfitr, 9± = 
sin 2 6\ , and e is the detector efficiency. The efficiency 
moments are determined using simulated events. The 
efficiency moments are fit to second-order even polyno- 
mials in cos#tr, the parameters of which are fixed in the 
subsequent likelihood fit to the cos 9 tI distribution. 

The measurement of R± is based on a combined un- 
binned maximum likelihood fit of the cos0 t r and tues 
distributions. The probability density function (pdf) for 
the ?71es distribution is given by the sum of ARGUS and 
Gaussian functions. The background shape is modeled 
by an even second-order polynomial in cos#tr- The pdf 
used for signal events is given by Eq. ^ The experimen- 
tal resolution of 9tr is not negligible and is accounted 
for by convolving the signal pdf with a double Gaussian. 
Also, the resolution of #tr has significant tails caused by 
mis-reconstructed events. The effect of these tails is ac- 
counted for by an additional term in the signal pdf. The 
parameterization of the 9t T resolution is determined from 
simulations. 

We categorize our events in three types: D* + D*~ — > 
(D a n+,D°7T-), (D a 7T+,D-TT°), and {D+ir°, D°^) be- 



cause events with a neutral slow pion and events with a 
charged slow pion have different background levels, de- 
tection efficiencies, and cos 6> tr resolutions. Thus, the pa- 
rameters determined in the likelihood fit are three signal 
fractions, the cos 9t r background shape parameter, three 
toes parameters (a and mean of the Gaussian, and k from 
the ARGUS function), and R±. The fit to the dataset 
yields a value of 

R± = 0.063 ± 0.055(stat) ± 0.009(syst). (3) 

Figure^shows the distribution of cos 9t r for events in the 
range toes > 5.27 GeV/c 2 . The value of a is fixed to zero 
in the fit, incurring a (negligible) systematic uncertainty. 
The largest systematic uncertainties arise from the pa- 
rameterization of the angular resolution (0.005) and the 
determination of the efficiency moments (0.005). 
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FIG. 1: Measured distribution of cos 6 tr and fit results. The 
data points are from the region mes > 5.27 GeV/c 2 and the 
solid line is the projection of the fit result in the same region. 
The dotted line represents the background component. 



In addition to the time-independent measurement of 
the CP-odd fraction, we perform a combined analysis of 
the cos 9tr distribution and the time dependence in order 
to determine the time-dependent CP asymmetry, using 
the sample of B° — > D* + D*~ events described previ- 
ously. We also use information from the other B meson 
in the event to tag its flavor as either a B° or B°. 

Although factorization models predict a small penguin 
contamination in the weak phase difference in Im(Ay) = 
— sin2/3 Q , a sizable penguin contribution cannot a pri- 
ori be excluded. Thus, the value of A/ = 7 7cPp^jj Q 
can be different for the three transversity amplitudes 
(/ =-L,0, ||) because of possible different penguin-to-tree 
ratios. This possibility is explicitly included in the pa- 
rameterization of the decay rates described here. 

The decay rate P+(P_) for a neutral B meson tagged 
as a B°(B°) is given by 

-|At|/-r B o r 

P ± (At) = — (G(l-lAP)T 



4r £ 

V[Ssm(Am d At) 



Ccos(ATO d A£)]},(4) 
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where At — t ICC — t tag is the difference between the proper 
decay time of the reconstructed B meson (B rcc ) and of 
the tagging B meson (B tag ), t b o is the B° lifetime, and 
Amd is the mass difference determined from the B°-B a 
oscillation frequency. The dilution factor, T> = 1 — 2ui, 
where Li is the average mistag fraction, describes the ef- 
fect of incorrect tags, and AT> accounts for possible dif- 
ferences in the mistag probabilities for B° and B°. The 
G, C and S coefficients are defined as 



G= §[(1-P ± )sin 2 



S = 



[(1 



l+IA+l 2 
f>, \ 2Im(X+) 
U ^> 1+|A+ 



2R ± cos 2 

,2 , 



(5) 



C = \ [(1 - i?±)^TF Sin2 tr + 2^^^ COS 2 tr ], 



^sin 2 fltr-2iix^ 



+ |Aj 



Because the two CP-even transversity amplitudes pro- 
duce the same distribution in cos (9 tr , we are only sensitive 
to A + , the appropriate average of Am and Xq: 



MA+) 
1 + IA+l 2 

1-IA+l 2 
1 + IA+l 2 



Im (A||) M 2 

WW W 



+ 



(Aol M 2 



l+|Ao 



MA,, 



Mi 



Ml 

l-|Ao ' 



I -MA \- "| T 1 + |A |2 M 



M 2 



M 2 



(6) 



If angular acceptance effects are not taken into account 
and the CP-odd fraction is allowed to float in the fit, then 
no bias is seen in the resulting value of A+ based on sim- 
ulations. Hence, a dedicated method to correct for detec- 
tor efficiency is not required. The value of R± obtained 
is therefore an effective value, which is not identical to 
the acceptance-corrected value from the time-integrated 
measurement. 

The time interval At is calculated from the measured 
separation Az between the decay vertex of the recon- 
structed B meson and the vertex of the flavor-tagging 
B meson along the collision axis. Events with a At un- 
certainty < 2.5 ps, and a measured \ At\ < 20 ps are ac- 
cepted. The mistag fractions and At resolution functions 
are determined from a sample of neutral B decays to fla- 
vor eigenstates, Bfl av , as in the sin2/3 measurement using 
charmonium decays |9j. Vertex reconstruction, the deter- 
mination of At, and the algorithms used for the determi- 
nation of the flavor of P tag are described in Refs. 0, 0] . 

We determine the parameters Im(A + ) and |A + | with 
a simultaneous unbinned maximum likelihood fit to the 
At distributions of the B rac and Bfl av tagged samples 
(Fig. El . The At distribution of the Bfl av sample evolves 
according to the known frequency for flavor oscillations 
in neutral B mesons. The observed magnitude of the CP 
asymmetry in the B lcc sample and the flavor oscillation 
in the Pfl av sample are reduced by the same factor T> due 
to flavor mistags. The At distributions for the B lec and 
-Bfl a v samples are both convolved with a common At res- 
olution function. The 9 tT angular resolution is accounted 
for in the same way as described previously. Events are 



assigned signal and background probabilities based on 
their mES values. Backgrounds are incorporated with an 
empirical description of their At distributions, contain- 
ing prompt (zero lifetime) and non-prompt components 
convolved with a separate resolution function |jj . 

A total of 38 parameters arc varied in the fit: the val- 
ues of Im(A+) and |A+| (2), the effective CP-odd frac- 
tion (1), the average mistag fraction w and the difference 
Aw between B° and B° mistags for each tagging cat- 
egory (8), parameters for the signal At resolution (9), 
and parameters for the background time dependence (7), 
At resolution (3), and mistag fractions (8). Because the 
CP-odd fraction is small, we have little sensitivity to the 
parameters |Aj_| and Im(Aj_). Therefore they are fixed 
to 1.0 and —0.741 respectively. These are the values 
expected if direct CP violation and contributions from 
penguin diagrams are neglected. The changes in the fit- 
ted values of Im(A+) and |A+| for different input values of 
Im(Aj_) (varied between —1.0 and 1.0) and |Aj_| (varied 
between 0.7 and 1.3) are taken into account as systematic 
uncertainties. The results obtained from the fit (Fig. EJ) 
are 

Im(A + ) = 0.05 ± 0.29(stat) ±0.10(syst) (7) 
|A+| = 0.75 ± 0.19 (stat) ± 0.02 (syst). (8) 

The dominant sources of systematic uncertainty come 
from the variation of the value of A^ (0.056 and 0.008, 
respectively, for Im(A+) and |A+|), and the level, compo- 
sition, and CP asymmetry of the background (0.078 and 
0.005). 

If the B —> D* + D*~ transition proceeds only through 
the b — > ccd tree amplitude, we expect that Im(A + ) = 
— sin 2/3 and |A+| = 1. To test this hypothesis, we fix 
Im(A + ) = -0.741 9] and |A+| = 1 and repeat the fit. 
The observed change in the likelihood corresponds to 2.5 
standard deviations (statistical uncertainty only). 

In summary, we have reported a measurement of the 
CP-odd fraction and measurements of time-dependent 
CP asymmetries for the decay B° — > D* + D*~. The 
measurement of R± supersedes the previous BABAR re- 
sult 01 > with a factor of three reduction in the statistical 
uncertainty, and indicates that B° — » D* + D*~ is mostly 
CP-even. The time-dependent asymmetries arc found to 
differ slightly from Standard Model predictions with pen- 
guin amplitudes ignored. 
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FIG. 2: From top to bottom: Number N B o (N S o) of candi- 
date events in the region mES > 5.27 GeV/c 2 with a B° (B°) 
tag, and the raw asymmetry (N B o — Ng )/(N B o + Ago), as 
functions of At. The solid curves represent the result of the 
combined fit to the full sample. The shaded regions represent 
the background contributions. 
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